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Thrust  Efficiency,  Energy  Efficiency  and  the  Role  of  VDF  in 
Hall  Thruster  Performance  Analysis*  (Preprint) 


C.  William  Larson1,  Daniel  L.  Brown2,  and  William  A.  Hargus,  Jr.3 

Air  Force  Research  Laboratory,  Edwards  AFB,  CA  93524-7680 

[Abstract]  A  rocket  power  efficiency  equation  was  written  to  explicitly  account  for  the 
effect  of  the  velocity  distribution  function  (VDF)  of  the  propellant  jet  on  the  conversion  of 
anode  electrical  energy  to  jet  kinetic  energy.  This  enabled  a  mathematically  rigorous 
distinction  to  be  made  between  thrust  efficiency  and  energy  efficiency.  In  this  approach 
anode  thrust  efficiency,  the  commonly  reported  figure  of  merit  for  Hall  thrusters,  is  the 
product  of  three  utilization  efficiencies:  (1)  propellant  utilization  efficiency,  (2)  voltage 
utilization  efficiency,  and  (3)  current  utilization  efficiency,  which  are  less  than  unity  under 
all  real  operating  conditions.  Unit  propellant  utilization  is  characterized  by  100  percent 
ionization  to  a  single  ionic  species  whose  thrust  vectors  are  all  directed  along  the  same  thrust 
axis.  Anode  voltage  utilization  efficiency  is  unity  when  ion  species  are  created  at  the  anode 
and  accelerated  through  the  entire  anode  potential.  Current  utilization  efficiency  is  the 
fraction  of  cathode  electron  flow  utilized  in  neutralization  of  accelerated  positive  ions.  It  can 
never  be  unity  because  a  portion  of  the  electron  flow  must  be  recycled  back  to  the  anode  to 
provide  energy  to  ionize  neutral  propellant.  The  architecture  of  the  efficiency  analysis  is 
such  that  energy  efficiency  becomes  naturally  expressed  as  a  product  of  voltage  and  current 
utilization  efficiencies,  and  is  rigorously  separated  from  propellant  utilization  efficiency. 
Thus,  thrust  efficiency  is  the  product  of  propellant  utilization  efficiency  and  energy 
efficiency.  The  methodology  is  applied  to  analysis  of  data  from  systematic  low  and  high 
power  (0.2  to  50  kW)  Hall  thruster  performance  studies  published  in  the  open  literature. 
The  cited  data  includes  measurement  of  thrust,  propellant  mass  flow  rate,  anode  voltage, 
and  anode  current  coupled  with  various  electrical  and  optical  diagnostics  that  provide 
information  about  the  VDF  and  thermodynamic  state  of  the  propellant  jet.  At  their 
optimum  operating  points  all  Hall  thrusters  appear  to  require  ~  20%  to  30%  of  the  total 
electron  flow  to  be  recycled  through  the  plasma  to  ionize  propellant  and  sustain  the  plasma 
processes.  Recycled  electrons  loose  their  energy  by  ionizing  neutrals  and  by  joule  heating.  It 
appears  that  anode  thrust  efficiencies  in  the  neighborhood  of  75%  that  have  been  achieved 
in  higher  power  thrusters  (>  ~  5  kW,  operating  with  xenon  and  producing  thrust  to  power 
ratio  ~  50  mN/kW  at  specific  impulse  ~  3000  seconds),  may  be  nearing  a  practical  upper 
limit  as  this  requires  a  geometric  mean  of  the  three  utilization  efficiencies  of  around  91%. 
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Nomenclature  (In  order  of  usage) 

rj thrust  —  experimental  thrust  efficiency. 

F  =  measured  thrust  vector. 

Va  =  measured  anode  voltage 

Ia  =  measured  anode  current 

Pin  =  input  power  to  anode. 

m  =  measured  propellant  mass  flow  rate 

<v>  =  mass  weighted  average  exit  velocity,  specific  impulse. 

g0  =  Earth’s  gravitational  constant  at  sea  level,  9.806  m/s2. 

Isp  =  specific  impulse,  independent  of  unit  system,  Isp  =  <v>/gG. 

uin  =  specific  energy  input  =  Pin/  m  =  Pin/[F/<v>]  =  <v>/[F/Pin], 

<v2>  =  mass  weighted  average  squared  exit  velocity. 

Pjet  =  power  of  jet  =  Vi  rii  <v2>. 

q energy  =  energy  efficiency  =  Pjct/Pin- 

<P  =  Phi  Factor  =  <v>2/<v2>,  propellant  utilization  efficiency. 

<Pvdf  =  Phi  factor  for  VDF  loss. 

®div  =  Phi  factor  for  divergence  loss. 

a  =  plume  momentum  divergence  half  angle. 

fo,fi,f2,f3  =  exit  mass  fractions  of  Xe°,Xe+1,Xe+2,Xe+3.  ti,  -  tj  +  t'2  +  f  =  1 . 

v0,v i ,v2,v3  =  exit  velocities  of  Xe°,Xe+1,Xe+2,Xe+3. 

r0,r2,r3  =  velocity  ratios,  v0/vi,  v2/vi,  v3/vi. 

fi  =  ion  mass  fraction  at  exit  =  f|+f2+f3. 

f*i,f  2,f*3  =  reduced  ion  mass  fractions  at  exit,  f*i+f* 2+f* 3=1,  f'l^f/f .  f* 2=f2/fi,  f* 3=f3/fi- 
a  =  quality  factor  of  VDF,  a  =  <Pvd/fi- 

3  =  Faraday  constant,  96485  coulombs/mol  of  charge. 

M  =  molecular  weight  of  propellant,  xenon  =  0.13 13  kg/mol,  krypton  =  0.08380  kg/mol. 

<v2>+  =  average  squared  velocity  of  ion  species  =  f  i  v  1 2+f* 2  v22+f* 3  v32. 

Q  =  average  charge  of  ion  species  =  f* i  +  2f*2  +  3f*3. 

SVi,  5V2,  5V3=  acceleration  voltages  of  Xe+1,Xe12,Xe+3. 

AV  =  average  acceleration  voltage  =  (f* i  5V  +2f* 2  5V2+3f* 35V3)/Q. 

(l-(3)  =  average  voltage  utilization  efficiency  =  AV/Va. 

(i  =  average  fractional  loss  of  acceleration  voltage  =  1-  AV/Va. 

X  =  Chi  Factor  =  f;Q,  charge  fraction,  moles  of  charge  per  mol  of  propellant. 

(1-r)  =  current  utilization  efficiency, 

r  =  fractional  loss  of  current,  electron  recycle  fraction. 

rmjn  =  minimum  fraction  of  current  required  to  ionize  propellant 

e  i,£2,£3  =  xenon  ionization  potentials. 

Ej  =  first  experimental  parameter  group  =  !/2(F/  m  f/CV.JJ/M )  =  Of  1  -(/)•/ 

E2  =  second  experimental  parameter  group  =  ( m  =  ( I  -r)/% 

I.  Introduction 

Flail  thrusters  enable  in-space  missions  that  require  thousands  of  hours  of  thrusting  with  tens  to  thousands  of  mN 
of  thrust.1  Numerous  studies  of  their  performance  have  appeared  in  open  literature  during  the  last  fifteen  years. 
Measurement2"16  of  performance  is  reported  in  terms  of  thrust  efficiency  (qthmst),  defined  in  terms  of  measured 
quantities,  F,  Pin  and  m  by: 


thrust  2  Pin  m 

In  Equation  (1)  the  quantity  F/P;n  is  the  thrust  to  power  ratio,  an  important  figure  of  merit  in  Hall  thruster 
performance.  The  second  important  figure  of  merit,  F/m  ,  is  the  propellant  jet  specific  impulse  (impulse  per  unit 
mass),  which  has  dimensions  of  velocity  in  SI  units.  Formally,  it  is  the  mass  weighted  average  velocity  vector  of  the 
jet  species  at  the  exit  surface  of  the  thruster,  denoted  <v>.  The  term  “specific  impulse”  as  commonly  used  in  the 
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rocket  propulsion  literature  is  defined  by  Isp  =  <v>/g0,  where  gQ  is  the  Earth’s  standard  sea-level  gravitational 
constant,  g0  =  9.806  m/s2.  Pm  is  the  anode  power,  given  by  P;n  =  VaIa,  where  Va  and  Ia  are  the  measured  quantities. 
Laboratory  measurements  are  made  by  setting  constant  Va  and  m ,  and  then  measuring  F  and  Ia  after  tuning  the 
magnetic  field  strength  and  topology  to  minimize  Ia  (or  maximize  r|  thrust)-  Specific  energy  investment  is  defined  by 
Ujn  =  Pin/  rii ,  which  may  also  be  written  uin  =  <v>/[F/Pin].  Thus,  a  standard  graph  of  figures  of  merit,  e.g.,  thrust  to 
power  ratio  versus  specific  impulse,  may  include  lines  of  constant  uln  that  pass  through  the  chart  origin  (F/Pin  = 
<v>/uin)  as  well  as  lines  of  constant  pthmst  that  may  be  produced  from  Eq.  (1).  Finally,  we  note  that  thrust  efficiency 
may  be  expressed  by  functions  explicit  in  uin:  r|thmst  =  !/2Uin(F/Pin)  ,  and  pthmst =  1/4<v>2/u;n. 

Hall  Thrusters  have  uin  values  around  100  to  1000  MJ/kg,  about  10  to  100  times  larger  than  chemical  propulsion 
systems.  Their  range  of  specific  impulse  may  vary  by  a  factor  of  ten,  from  ~  3  to  ~  30  times  larger  than  found  in 
chemical  propulsion.  Thrust  efficiency  varies  drastically  between  low  and  high  power  Hall  thrusters:  r|thmst  <  0.5 
for  Pin  <  ~  600  W,  to  q thmst  ~  0.7  for  Pin  >  ~  5  kW.  High  power  thrusters  appear  to  be  near  their  practical  upper  limit 
of  thrust  efficiency.  One  major  determinant  of  thrust  efficiency,  the  “propellant  utilization  efficiency”,  is  defined  in 
at  least  three  different  ways  in  the  literature,  i.e.,  most  commonly  as  the  ionization  fraction,  but  also  as  the  fraction 
of  momentum  carried  by  ions1,  or  the  beam  current  charge  fraction2,  or  the  quantity  obtained  by  integrating  Faraday 
probe  measurements  of  beam  current  around  the  jet  axis  at  some  distance  from  the  exit.  In  this  paper  we  develop  a 
fundamentally  rigorous  and  simple  framework  that  reconciles  various  usages  of  propellant  utilization  and  other 
utilization  efficiencies  to  enable  consistent  analysis  of  various  aspects  of  performance,  and,  in  particular,  to  enable 
comparison  of  performances  of  low  power  Hall  thrusters  to  high  power  Hall  thrusters. 


II.  Velocity  Distribution  Function 

When  the  velocity  distribution  function  (VDF)  of  the  propellant  jet  over  the  exit  surface  (exhaust)  of  the  thruster 
is  measured,  an  experimental  energy  efficiency  (^energy  =  Pje/Pin)  may  also  be  extracted  from  the  measured 
quantities.  Since  Pjet  =  !4  m  <v2>,  where  <v2>  is  the  jet’s  mass  weighted  average  squared  velocity ,  an  exact  power 
efficiency  equation  may  be  written  that  applies  to  any  jet,  whether  continuous  or  pulsed  viz., 


1  F  Pjet 

n  = - <v>=  0^1  =  011 

'thrust  2  Pin  P:„ 


energy  ’ 


(2) 


where 


<D  = 


<  v> 


<v2> 


(3) 


Equation.  (2)  is  an  exact  mathematical  representation,  equivalent  to  <v2>qthrust  =  <v>2qenergy.  It  is  based  on  the 
First  Law  of  Thermodynamics  and  Newton  s  Second  Law,  Pjet  =  !4  m  <v2>,  and  -F  =  d(mv)/dt  =  m  <v>,  where  <v> 
is  the  specific  impulse  of  the  propulsive  device,  i.e.,  impulse  per  unit  mass  or  jet  velocity.  Knowledge  of  the  VDF 
(i.e.,  <P)  is  necessary  and  sufficient  for  calculation  of  the  energy  efficiency.  For  simplicity  in  the  current  discussion, 
the  Phi  Factor  may  be  factored  into  two  parts,  one  to  represent  loss  from  plume  divergence,  and  one  to  represent  loss 
from  non-uniformity  of  the  VDF  that  occurs  when  the  jet  is  composed  of  species  with  widely  differing  velocity. 

(I)  =  ®div®vdf  •  (4) 


For  an  isotropic  distribution  of  divergences,  with  maximum  momentum  divergence  half  angle  a,  <3>div  may  be 
expressed  by17 


O 


div 


1  +  cos  a 


(5) 


The  quantity  within  parenthesis  is  the  familiar  factor  by  which  thrust  of  an  ideal  1-dimensional  rocket  jet  is  reduced 
due  to  its  2-dimensional  divergence,  viz.,  F  =  m  <v>[(l+cos  a)/2] 
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III.  Phi  Factor  for  Simple  Velocity  Distribution  Functions 

The  <D  factor  is  mathematically  limited  to  0  <  <5  <  1.  For  a  mono  energetic  delta  function  distribution  of  velocities 
<Dvdf=  1.  For  a  Maxwellian  distribution,  M(v)  =  exp(-'/2  mv2/kT),  ®vdf  =  8/(3  ji).  For  a  Gaussian  distribution,  G(v)  = 
l/(27ia)1/2exp-[(v-<v>)2/(2c2)],  ®vdf  =  1/[1  +  (cr/<v>)2],18  where  a  is  the  standard  deviation.  Chemical  rockets 
produce  shifted  Maxwellian  distributions  that  are  similar  to  Gaussians  with  <v>  ~  3000  m/s  and  a  ~  <  10%  of  <v>, 
for  which  <Dvdf  >  0.99.  Thus,  the  <Dvdf  factor  does  not  often  appear  in  textbook  equations  where  scalar  energy 
quantities  are  compared  to  squared  vector  momentum  quantities. 

IV.  Trimodal  Ion  Distribution  with  Neutral  Xenon 

To  show  how  the  Phi  Factor  depends  on  the  properties  of  a  jet  with  non-uniform  velocity  distribution,  consider  the 
jet  that  contains  Xe°,  Xe+1,  Xe+2,  and  Xe+J  with  mass  fractions  of  f0,  fh  f2,  and  f3,  respectively,  so  that  f0+fi+f2+f3=l. 
To  further  simplify,  we  assume  the  neutral  and  each  ion  to  have  a  delta  function  distribution  of  velocities,  specified 
by  Vo,  Vi,  v2,  and  v3.  This  assumption  will  not  significantly  affect  conclusions  that  would  result  from  a  more 
complicated  model,  such  as  a  series  of  Gaussians  for  example.  In  this  case  the  VDF  Phi  Factor  is: 

(V0f0+V1f,  +  V2f2  +  V3f3  f  (6) 

vo2fo  +  Vj2  f  j  +  V22  f2  +  V2  f3 

For  convenience  and  for  compact  notation,  the  ion  mass  fraction  is  defined  by  f;  =  f3  +  f2  +  f3,  so  that  f(l  +  f;  =  1, 
and  reduced  ion  mass  fractions  are  defined  by  fj  =  fj/fj,  f2  =  f2/f;,  and  f*3  =  f3/f;,  so  that  f ,  +  f*2  +  f*3  =  1 .  Ratios  of 
velocities  of  Xe°,  Xe+2,  and  Xe+3  to  the  velocity  of  Xe+1  are  defined  by  r0  =  v0/vi,  r2  =  v2/V|,  r3  =  v3/vi.  This  enables 
<Svdf  to  be  written  as  a  function  explicit  in  any  two  of  the  three  reduced  ion  fraction  variables  (f/  ,  f2*,  and  f3  ).  In 
terms  of  fi  and  f2  the  result  is: 


ft) 

fl-fi 

l  f, 

+  1-3+  f*  (l  -  r3)  +  f*(r2-r3) 

2 

ro 

— 1 

f,  J 

+  r32  +  f*  (l _  r3  )  +  f*(r22  -  r32 ) 

(7) 


First  we  use  Eq.  (7)  with  f2*  =  0  and  r3  =  0  to  illustrate  the  dependence  of  <Pvdf  on  f  and  r0  for  a  two-component 
jet  composed  of  slow  moving  neutrals  and  fast  moving,  singly  charged  ions.1  Figure  1  shows  <Pvdf  as  a  function  of 
Vi/vo  for  this  bimodal  distribution,  with  values  of  f  from  0.1  to  1  in  increments  of  0.1.  The  Figure  shows  that  <Pvdf 
approaches  f,  as  Vi/v0  approaches  infinity;  when  f  is  larger,  as  Vi/v0  increases,  <D  approaches  f,  faster.  Based  on 
experimental  measurements,12  v0  ~  200  m/s,  and  Vi  ~  20,000  m/s,  so  r0  ~  0.01  is  appropriate  and  Ov(ir  ~  fi- 


Fig.  1.  Dependence  of  d>vdf  on  velocity  ratio  Vi/v0  =  l/r0  for  bimodal  velocity  distribution 
with  f|  values  from  0.1  to  1.0  in  increments  of  0.1. 
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Additional  thrust  inefficiency  occurs  when  the  jet  is  composed  of  three  ion  species.  This  effect  may  be  easily 
analyzed  for  an  idealization  where  all  three  species  are  created  in  the  same  zone  whose  length  is  small  compared  to 
the  acceleration  length,  in  this  approximation  ion  kinetic  energies  are  proportional  to  their  charge  so  that  r2  =  V2/V1  = 
a/2,  and  r3  =  v3/vi  =  a/3.  Additional  approximation  is  implicit  here  because  ionization  potentials  are  ignored  here, 
which  is  valid  so  long  as  they  are  a  small  fraction  of  the  acceleration  (anode)  voltage.  Figure  2a  shows  <Dvdf  as  a 
function  of  f  jwith  these  velocity  ratios  for  all  f*1:f*2:f*3  distributions  with  f;  values  from  0.1  to  1  in  increments  of 
0.1.  Figure  2b  shows  detail  of  the  f;  =  0.9  manifold.  Experimental  findings  from  RPA  and  mass  spectrometric 
probes  of  various  thrusters  are  also  plotted  in  Fig.  2b.  Ekholm  ,  et  al 1 5,  reported  [fpfpf*^  [78:15:7],  for  a 
relatively  highly  charged  jet  (average  charge,  Q  =  f*]+2f*2+3f*3  =  1.29)  of  a  Busek  Hall  Thruster  (BHT  600) 
running  under  optimum  conditions  with  Pin  =  600  W,  Va  =  300V,  m  =  2.5  mg/s,  r|thrast  =  0.555. 14  Hofer,  et  al.4,5 
studied  several  NASA  Hall  Thrusters.  For  the  173Mv2  thruster,  they  reported  variation  of  [f  pf* 2:f  3]  from  [96:3:1], 
Q  =  1.05,  to  [81:16:3],  Q  =  1.22  as  Va  increased  from  300  to  900  volts,  and  specific  impulse  increased  from  ~  10  to 
30  km/s.  Gulzcinski,  et  al.10  made  mass  spectrometric  measurements  in  the  near  and  far  field  of  plumes  of  the  P5 
Hall  thruster  and  reported  large  decreases  in  average  charge  as  distance  from  the  exit  increased;  their  measurement 
for  the  distribution  with  largest  Q  is  indicated  in  Fig.  2b.  The  departure  of  <J>vd(  from  f  at  f;  =  0.9  is  less  than  1%  to 
2%  for  Hofer  distributions  and  ~  3.5%  for  the  Ekholm  and  Gulczinski  distributions. 


Fig.  2.  (a)  VDF  Phi  Factor  for  jet  composed  of  neutrals  and  three  ions.  Ten  manifolds  of  f  l5  f  2,  f  3 
distributions  are  shown  at  f;  values  from  0. 1  to  1  at  f;  increments  of  0. 1 . 

(b)  Detail  of  <J>vdf  for  f;  =  0.9  manifold  showing  lines  of  constant  f  2  (blue  solid  lines)and  constant  f  3 
(black  dashed  lines)  from  0  to  1  and  lines  of  constant  Q  (green  dotted  lines)  from  1.1  to  1.6.  The  locations 
of  points  for  the  Ekholm  distribution,  several  Hofer  distributions  (where  <J>vdf  decreases  with  increasing  Q 
and  Va),  and  the  Gulczinski  distribution  are  also  shown. 


Figure  3  shows  the  deviation  of  <Pvdf  from  the  <Pvdf  =  f,  approximation  as  a  graph  of  [  ®vdr  -  fi]  vs  t).  Eq.  (7)  may 
be  used  to  show  that  the  envelope  of  all  deviations  is  bounded  by  the  hypothetical  distributions  [f  pf  2:f  3]  =  100:0:0 
(Q  =  1)  and  [60:0:40]  (Q  =  1.8).  The  Figure  also  shows  the  [<3>vdf  — fj]  vs  f,  lines  for  Hofer  and  Ekholm  distributions. 
It  is  apparent  that  [<Pvdf  -  fi]  may  be  expressed  by  a  linear  function  of  f;  over  a  large  range  of  fj  and  that  the  effect  of 
VDF  non-uniformity  may  be  approximated  with  a  simple  relationship:  d>vdf  ~  at).  Table  1  summarizes  details  for 
these  distributions  and  establishes  a  highly  correlated  linear  relation  between  Q  and  a,  i.e.,  Q  =  9. 6-8. 6a. 


Table  1.  Summary  of  ion  distributions  in  three  Hall  thrusters:  Q  =  9. 6-8. 6a. 


distribution 

fiHH* 

Q 

a  =  <Pvdf/fi  from  Eq. 
(7)  with  f  =  0.9 

limit 

100:00:00 

1.00 

1.000 

Hofer  173Mv2 

96:03:01 

1.05 

0.993 

Hofer  173Mv2 

93:06:01 

1.08 

0.989 

Hofer  173Mv2 

87:12:01 

1.14 

0.982 

Hofer  173Mv2 

81:16:03 

1.22 

0.973 

Ekholm  BHT  600 

78:15:07 

1.29 

0.963 

Gulczinski  P5 

70:23:07 

1.37 

0.961 

5 

American  Institute  of  Aeronautics  and  Astronautics 


0.0  0.2  0.4  0.6  0.8  1.0 

fi 

Fig.  3.  [®vdf—  fj]  vs  fj.  Departure  of  Phi  Factor  from  <J>vdf  =  fj  approximation  as  a 
function  of  fj  for  hypothetical  distributions  and  Ekholm,  and  Flofer  distributions. 


V.  Propellant,  Voltage,  and  Current  Utilization  Efficiencies 

Equation  (1),  may  be  written  in  a  form  that  reveals  thrust  efficiency  to  be  a  product  of  three  separate  efficiency 
terms,  viz., 

,  2 


1  2 

<v2>  2<V~>+ 


thrust 


2  7 

<  V  >  <\z  > 


+ 


— V 

M  a 


m  5 

TaM 


(8) 


The  first  term  in  Eq.  (8)  is  the  O-factor,  which  is  defined  by  Eqs.  (2-3).  As  used  in  this  paper,  the  O-factor  has 
attributes  of  a  good  “propellant  utilization  efficiency.”  As  Eq.  (3)  shows,  it  is  always  less  than  unity  and  it 
approaches  unity  as  propellant  utilization  approaches  perfection.  It  captures  the  loss  of  propellant  utilization  in  three 
ways:  (1)  by  quantifying  loss  of  utilization  due  to  incomplete  ionization,  (2)  by  quantifying  loss  of  thrust  due  to 
distribution  of  energy  among  several  ion  species  as  in  Eqs.  (6-7),  and  (3)  by  quantifying  loss  of  thrust  due  to 
divergence  of  the  propellant  trajectory  from  the  thrust  axis  as  in  Eqs.  (4-5).  The  <D-factor  is  unity  for  a  perfect, 
mono  energetic  jet  with  one  ion  species,  no  neutral  atoms,  and  no  divergence.  Also,  as  shown  by  Fig.  3,  Ovdf  is  only 
slightly  less  than  the  ionization  fraction.  Several  usages  of  the  term  “propellant  utilization”  have  appeared  in  the 
literature,  ranging  from  the  review  of  Zhurin  et  al.,1  where  it  is  the  fraction  of  momentum  earned  by  ions,  to  the 
paper  by  Raitses,  et  al.2  where  it  is  the  ratio  of  beam  current  to  mass  flow,  to  frequent  usages  where  it  is  simply 
referred  to  as  the  ionization  fraction. 

The  second  term  of  Eq.  (8)  may  be  transformed  into  an  expression  that  contains  explicit  voltage  utilization 
efficiencies  for  each  ion  species,  5V,/Va,  SVj/Va,  5Vj/Va,  such  as  may  be  directly  measured  with  an  ExB  probe. 
Also,  an  expression  may  be  derived  that  defines  an  average  voltage  utilization  efficiency,  AV/Va,  such  as  may  be 
directly  measured  with  a  retarding  potential  analyzer,  RPA.  First,  one  may  obtain  the  leading  ratio  of  the  second 
term  of  Eq.  (8),  <v2>/<v2>+,  as  follows: 


<V2>  _  fpVo2  +  f lvf  +f~2v22  +^3V32  f 
<V2>+  ffvf  +f2*V22  +f3*v32  ' 


1  + 


f*  +f2*r22  +f3*r32 


fi 


(9) 


It  may  be  shown  that  the  bracketed  quantity  of  Eq.  (9)  is  very  close  to  unity  under  all  reasonable  conditions  of 
operation  of  a  Flail  thruster  because  the  term  in  numerator  of  the  quotient  is  very  small,  i.e.,  r02(l-fi)/fi  ~  lxlO"5, 
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compared  to  the  denominator.  Thus,  the  approximation,  <v2>/<v2>+  ==  t,  is  accurate  to  better  than  ten  parts  per 
million. 

The  kinetic  energy  of  an  ion  species  is  proportional  to  the  product  of  its  charge,  z,  and  acceleration  voltage,  5V, 
e.g.,  ‘/2V2  =  zoVT/jlt,  so  the  average  specific  ion  energy  may  be  written  for  a  propellant  jet  composed  of  three  ion 
species: 


-  <  v2  >  = — (  f *5V|  +  2  f  *  5V2  +  3f*  6V3 ) , 
2  +  M  '  1  22  3  i  ’ 


(10) 


where  5V[,  8V2,  and  5V3,  are  the  voltages  through  which  the  three  ions  are  accelerated.  Using  Eq.  (9)  and  (10),  one 
may  write  the  second  term  of  Eq.  (8)  as  follows: 


1  2 

9  —  <  v  > 

<v/>  2  + 


<  v2  >, 


q 

—  V 
M  a 


=  f: 


8Vl 

Va 


+  2f, 


§V2 

Va 


+  3U 


5V3 

Va 


(11) 


An  average  voltage  utilization  efficiency,  AV/Va  =  ( 1  -(i).  may  be  defined  so  as  to  make  use  of  measurement  of 
kinetic  energy  (voltage)  per  unit  charge  with  an  RPA.  Thus, 

=  f  ^(^+2^+3^)=(l-P)fjQ  (12) 

V  a 

In  Eq.  (12)  the  average  charge  is  defined  by  Q  =  f|*  +  2tV  +  3tY.  Average  acceleration  voltage  is  defined  by  AV  = 
[  f|  5Vi  +  t'2  5V2  +  f3  5V3]/Q.  The  ft  quantity  is  the  average  fractional  acceleration  voltage  loss,  so  that  ( 1  -[:’>)  = 
AV/Va  is  the  average  voltage  utilization  efficiency. 

Another  useful  property  of  the  jet  that  appears  in  Eq.  (12)  is  denoted  the  Chi-factor,  x  =  f\Q,  which  is  similar  to 
the  O-factor  because  they  are  both  proportional  to  f.  Although  the  '/-factor  contains  information  about  ionization 
fraction,  it  is  not  be  suitable  as  a  utilization  efficiency  because  it  may  have  a  value  greater  than  unity,  e.g.,  in  highly 
ionized  plasmas  with  highly  charged  ions.  The  quantity  [/  m  J/M]  is  the  total  positive  ion  current  exiting  the 
thruster,  frequently  called  the  beam  current,  which  is  the  quantity  obtained  in  a  good  Faraday  probe  measurement.  7 
may  also  be  thought  of  as  the  positive  charge  fraction,  or  as  the  moles  of  positive  charge  flow  per  mole  of  propellant 
flow. 

The  third  term  of  Eq.  (8)  may  be  transformed  into  the  quotient  of  a  current  utilization  efficiency,  ( 1  -r),  and  the  7  - 
factor,  viz., 


SVi  *8^+3f.5V3 

1  \T 


Va 


Va 


va 


m  5  _  (1-r) 
Ia  M  X 


(13) 


by  transforming  the  equation  of  mass  and  charge  conservation,  viz., 

(1  -  r)Ia  =  m  -^-x  =  m  f,  (f*  +  2f2*  +  3f3*).  (14) 

In  Eqs.  (13)  and  (14),  r  is  the  “electron  recycle  fraction”,  the  fraction  of  electrons  that  recycle  from  the  cathode, 
through  the  acceleration  channel,  and  to  the  anode.  Recycle  electrons  ionize  neutral  propellant  and  sustain 
dissipative  plasma  processes.  Energy  not  used  for  ionization  is  lost  by  joule  heating;  the  total  power  loss  is  rIaVa. 
The  rla  quantity  may  be  thought  of  as  the  Hall  current  required  to  sustain  plasma  ionization.  The  quantity  (l-r)Ia 
represents  the  beam  current,  which  is  the  electron  current  utilized  in  neutralization  of  the  positive  charge  flow 
leaving  the  thruster.  The  minimum  electron  recycle  fraction,  rmin,  is  the  fraction  of  the  anode  current,  la,  that  carries 
the  amount  of  power  necessary  to  produce  a  given  ionization  fraction  of  specified  ion  distribution,  which  is  given 

by: 
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m 


(15) 


.  fi(eifl*+  £2f2*+e3f3*) 

^  a 

where  £1,  £2,  and  £3  are  the  first,  second  and  third  ionization  potentials  of  xenon.  Eq.  (15)  may  also  be  written  as  rmjn 
=  uion/Ujn,  where  uion  =  f;(£ifi  +  £2tV  +  S3f3  )  is  the  specific  energy  required  to  produce  ion  fraction  f,  with  specified 
species  distribution  [fi  :f2  :f3  ]. 

Substitution  of  the  mass,  energy  and  charge  conservation  principles  embodied  in  Eqs.  (9)-(14)  into  Eq.  (8)  allows 
thrust  efficiency  to  be  written  as  a  product  of  three  utilization  efficiencies,  viz., 

Vrast=°(1-P^  *^-^=0(1-p)(1-r),  (16) 

and  the  energy  efficiency  to  be  written  as  a  product  of  voltage  and  current  utilization  efficiencies,  independent  of, 
and  rigorously  separable  from,  propellant  utilization  efficiency,  viz., 

n  =  (1  —  P)(  1  —  r) .  (17) 

energy 

The  recent  paper  by  Hofer  and  Gallimore11  presents  a  comprehensive  study  of  the  performance  of  the  NASA 
173Mv2  Hall  thruster  and  sets  forth  an  approach  for  efficiency  analysis  that  includes  definitions  of  voltage 
utilization  efficiency  ( 1  -fi)  and  current  utilization  efficiency  (1-r)  that  are  identical  to  the  definitions  used  here. 
However,  in  their  system,  anode  thrust  efficiency,  defined  as  in  Eq.  (1)  of  this  paper,  is  the  product  of  four 
utilization  efficiencies,  (l-(3),  (1-r),  and  two  other  utilization  efficiencies,  termed  mass  utilization  efficiency,  r|m,  and 
charge  utilization  efficiency,  r|q.  The  product  of  their  mass  and  charge  utilization  efficiency,  qqqm,  is  almost  exactly 
equivalent  to  our  <t>vdr  =  af;  as  defined  by  Eq.  (7),  and  as  approximated  by  fl>vdf =  af,  where  f,  represents  the  quantity’ 
of  propellant  utilization  and  a  represents  the  quality  of  propellant  utilization.  The  Hofer-Gallimore  qrn  is  identical  to 
f.  and  their  qq  is  nearly  identical  to  our  a,  except  that  their  qq  does  not  include  the  terms  involving  neutral  xenon, 
i.e.,  terms  involving  r0  in  Eq.(7).  But,  this  difference  makes  their  qq  ~  0.2%  smaller  than  our  a  when  f;  >  O.9.. 
Further,  the  difference  between  a  and  r|q  disappears  entirely  when  the  propellant  ionization  fraction  reaches  unity. 

The  Hofer  Gallimore  expression  for  anode  thrust  efficiency  does  not  include  a  term  for  loss  of  thrust  due  to 
momentum  divergence  of  the  propellant  jet.  Eq.  (5)  shows  that  the  divergence  loss  factor  may  be  significant.  When 
the  momentum  divergence  half  angle  increases  from  a  =  20-to  30-degrees,  <l>d,v  decreases  from  0.94  to  0.87.  Thus, 
the  ionization  fraction  (fj)  achieved  by  the  173Mv2  that  was  reported  as  0.87  will  be  significantly  larger  when 
divergence  loss  is  accounted  for,  i.e.,  fj  =  0.93  if  a  =  20-degrees,  or  f  =  1.0  if  a  =  30-degrees.  Also,  since  the 
ionization  fraction  must  be  less  than  unity,  a  must  be  less  than  30-degrees.  Since  average  ion  charge  was  measured 
by  ExB  probe  (Q  =  1.2)  and  %  =  1.2  is  obtained  from  total  beam  current  measurement,  the  value  of  f;  =  0.95  may  be 
obtained. 


VI.  Extraction  of  Utilization  Efficiencies  from  Experimental  Quantities 

Two  quantities,  Ei  and  E2,  that  mix  the  utilization  efficiencies  with  the  /-factor  may  be  obtained  directly  from 
groups  of  experimentally  measured  quantities. 


and 


1 

E,  =  0(  1  ~  p)X  =  OdjvaQ(  1  -  p)f,2  =  — 

V  — 

a  M 


F 

rii 


(1-r)  (1-r)  m  5 


(18) 


(19) 


Thrust  efficiency  then  may  then  be  expressed  simply:  qthmst  =  E]E2.  Thus,  graphs  of  Ej  vs  E2  with  lines  of  constant 
0 thrust  may  be  prepared  to  enhance  insightful  analysis.  For  example,  as  f  increases  while  other  parameters  [®d;v,  a, 
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Q,  r,  and  [S  |  are  constant,  Ej  increases  proportional  to  the  square  of  f  and  E2  decreases  proportional  to  the  reciprocal 
of  fj. 

Experimental  studies  that  include  measurement  of  beam  current  [beam  current  =  ( 1  -r)Ia]  with  use  of  Faraday 
probe  scans  of  the  plume  around  the  thrust  axis  and  at  some  distance  away  from  the  thruster  exit  enable  extraction  of 
additional  performance  data  because  %  may  be  obtained  by  Eq.  (14).  Then  current  utilization  efficiency,  (1-r)  may 
be  obtained  directly  from  Eq.  (19),  and  the  product  of  propellant  utilization  efficiency  and  voltage  utilization 
efficiency,  ®(1-P),  may  be  obtained  directly  from  Eq.  (18).  Further  deconvolution  may  be  accomplished  by  using 
the  previously  discussed  relationship  between  Q  and  a  (Q  =  9.6  -  8.6a),  which  enables  setting  limits  on  permissible 
values  of  fi,  Q,  and  the  geometric  mean  of  the  momentum  divergence  phi  factor  and  voltage  utilization  efficiency, 
[a>div(l-p)]1/2.  Since  ti  =  %/Q  and 


>W®d,v(‘-P)  = 

a  graph  of  f  vs  [  'Ifni  I  -[>)  | ' 2  may  be  produced  that  shows  permissible  solutions  that  depend  on  Ej  and  Q. 

Direct  measurement  of  voltage  utilization  efficiency,  (1-  P),  may  be  accomplished  with  use  of  a  retarding 
potential  analyzer,  and  performance  may  be  further  specified  by  an  expression  of  f;  as  a  function  of  4>div.  Finally,  if 
the  VDF  is  determined  with  an  ExB  probe  or  by  LIF,  a  complete  set  of  performance  parameters,  Q,  t’.  cpdiv  and  (1-P) 
may  be  deduced. 


8.6 


QA9.6-Q 


(20) 
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VII.  Applications  to  Hall  Thrusters  -  Cylindrical  Hall  Thruster 

An  extensive  literature  on  Hall  thruster  performance  exists.  Here  we  analyze  several  sets  of  notable  data2"16  to 
illustrate  the  distinctions  between  thrust  and  energy  efficiency  and  the  expression  of  thrust  efficiency  as  a  product  of 
the  three  utilization  efficiencies.  First,  we  examine  the  performance  of  a  Cylindrical  Hall  Thruster  (CHT)  reported 
by  Raitses,  et  al.2.  Five  CHTs  with  various  geometries  were  operated  at  powers  up  to  ~  600  W.  Systematic 
measurement  of  F,  m ,  Ia,  Va,  and  7  (with  Faraday  probe)  were  reported.  Results  from  all  the  geometries  were 
transcribed  from  the  graphical  presentations  in  the  published  paper  and  tabulated  in  Table  1  of  the  Appendix  to  this 
paper,  along  with  calculated  values  of  r|,hmst,  F/Pin,  <v>,  <I>(  1  -P)y,  and  (l-r)/%.  Also,  since  7  was  measured,  values 
of  the  current  utilization  efficiency  (1-r)  and  the  product  of  propellant  and  voltage  utilization  efficiency,  0(1 -P), 
could  be  extracted.  Other  relationships  between  Odlv,  f),  Q  and  their  limits  are  also  derived  from  the  data.. 

Here,  we  analyze  the  optimum  CHT  geometry,  which  had  a  30-mm  acceleration  channel  with  10-cm2  cross 
section.  Figure  4  shows  two  performance  runs  at  each  of  three  anode  voltages,  200  V,  250  V,  and  300  V.  Each  run 
consisted  of  a  thrust  measurement  at  each  of  six  flow  rates.  After  setting  Va  and  m ,  the  magnetic  fields  were 
adjusted  to  minimize  Ia.  Then  Ia  and  F  were  recorded.  At  the  same  time  the  total  positive  ion  beam  current  was 
measured  with  a  Faraday  probe,  and  this  enabled  determination  of  7.  The  7-factor  is  moles  of  positive  charge  per 
mole  of  propellant,  and  may  be  thought  of  as  a  molar  charge  fraction  as  defined  in  Eq.  (14). 


Fig.  4.  Performance  of  the  optimum  CHT  geometry  (30  mm  acceleration  channel  with  10  cm2  cross  section)  calculated  from  the 
data  reported  by  Raitses,  et  al.2.  Two  series  of  data  were  acquired  at  each  voltage:  Va  =  200  V  (squares),  250  V  (triangles),  and 
300  V  (circles).  Data  points  at  6  mass  flow  rates,  (0.88,  1.04,  1.29,  1.53,  1.84,  2.16  mg/s)  show  increasing  efficiency  with 
increasing  mass  flow  rate.  Lines  of  constant  thrust  efficiency  (rithmst =  0.2,  0.3,  0.4,  and  0.5)  and  lines  of  constant  specific  input 
energy  (uin  =  170-,  220-,  and  270-J/mg)  are  shown.  Note  that  uln  =  Va(/7/dt)y/(l-r).  At  300  V,  the  figure  shows  that  F/P  vs  <v> 
approximately  tracks  the  constant  Uj„  =  270  J/mg  line;  thus,  at  300  V,  when  flowrate  increases,  %/(l-r)  remains  constant. 
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Figure  5  shows  the  E2  vs  Ei  graph  and  the  '/-factor  (f;Q)  vs  Ia  graph  as  obtained  from  experimental  data. 

The  left  graph  shows  Ej  vs  E2  with  lines  of  constant  r|thrust,  which  is  based  only  on  the  F,  m ,  Ia,  Va  data.  The 
abscissa,  Ei,  is  proportional  to  the  square  of  the  ion  mass  fraction,  ff\  i.e.,  Ei  =  ®(  1  -fi)/  =  rl»dlvaf,(  1  -[i)f,Q-  The 
ordinate,  E2,  is  proportional  to  the  reciprocal  of  f;,  i.e.,  E2  =  ( I  -r)//  =  (l-r)/f;Q.  Thus,  the  pathway  of  increasing  t’  at 
constant  r,  fl.  a,  Q,  and  4>div  is  downward  to  the  right,  as  shown  in  the  left  graph.  The  Figure  shows  that  the  lower 
flow  rate  data  at  Va  =  200  V  and  300  V  follow  lines  of  increasing  f;  and  level  off,  coincident  with  Va  =  300  V,  at  (1- 
r)/y  ~  0.8.  The  /  vs  discharge  current  graph  on  the  right  side  of  Fig.  5  is  obtained  from  the  beam  current  as 
measured  with  a  Faraday  probe,  i.e.,  beam  current  =  (l-r)Ia  =  m  y(.7/Jt ). 


Ei  =  a>(1-p)x  =  <PdivaQ(l-P)fi2  Discharge  Current  (la,  A) 


Fig  5.  Performance  of  the  optimum  geometry  of  the  Cylindrical  Hall  Thruster  (CHT,  L  =  30  mm,  A  =  10  cm')  calculated 
from  data  of  Raitses,  Ashkenazy,  and  Guelman. 2.  Graph  of  Chi  factor  vs  discharge  current  was  obtained  from 
measurement  of  beam  current  by  integration  of  Faraday  probe  data. 


As  shown  in  Fig.  6,  use  of  the  /-factor  with  Eqs.  (18)  and  (19)  enables  construction  the  graph  of  (1-r)  vs  O(l-P). 
At  the  optimum  operating  point  the  measured  quantities  were  Va  =  300  V,  m  =2.16  mg/s,  Ia  =  2.02  A,  F  =  36.0 
mN,  and  /  =  0.881,  which  enabled  calculation  of  qthn,st  =  0.496,  E,  =  0.631,  E2  =  0.786,  (1-r)  =  0.692,  and  d>(l-P)  = 
0.716,  see  appendix.  Making  use  of  the  relationship  between  Q  and  a  (Table  1)  and  the  definitions  d>vdf=  afj  and  /  = 
fiQ,  and  Eq.  (18),  one  may  derive  a  relation  to  enable  consistent  sets  of  values  of  Q,  (1-P),  Odlv,  and  f  to  be 
determined. 

Equation  (20)  was  used  to  produce  Figure  7,  which  shows  a  plot  of  the  geometric  mean  of  the  divergence  part  of 
propellant  utilization  efficiency,  ®div,  and  the  voltage  utilization  efficiency  (1-P)  as  a  function  of  f,  for  each  of  the 
six  flow  rates  at  each  of  the  three  voltages  and  over  a  range  of  average  ion  charge,  Q  from  1.0  to  ~  1.35..  The  figure 
shows  the  restrictions  on  the  parameters  that  are  imposed  by  the  data  and  by  the  requirements  that  Q  >  1,  and 
[4>div(l-p)]‘/2  <  1.  Several  conclusions  may  be  drawn  from  Fig.  7.  For  example,  at  anode  voltage  of  300V  and 
flowrate  of  2.16  mg/s,  the  value  of  f  must  lie  between  ~  0.88  (where  Q  =  1,  [®div(l-P)]1/2  =~  0.9)  and  ~  0.73  (where 
Q  =  1.2,  [<I>div(l-P)]1 12  ~  1).  If  Q  remains  constant  at  Q  =  1.1  during  the  increase  of  flow  rate  from  0.88  to  2.16  mg/s, 
1’  must  increase  from  ~  0.60  to  0.80  and  [®div(l-P)]' 2  must  increase  from  0.86  to  0.95.  If  [(I)dlv(l-P)|12  remains 
constant  at  0.90  during  the  increase  of  the  flow,  fj  must  increase  from  0.56  to  0.88  and  Q  must  decrease  from  ~  1.18 
to  1.0.  Also,  at  the  highest  flow,  a  value  of  Q  =  1.1  is  consistent  with  f;  =  0.8,  (1-P)  =  .95,  and  4>div  =  0.95  (or  a  = 
18-degrees).  The  extreme  values  of  Q  =  1,  1)  =  0.88,  (1-P)max  =  0.965  that  apply  to  perfect  ionization  (i.e.,  f;  =  0.88, 
fi  =  1,  single  Xe'ion)  and  perfect  voltage  utilization  are  consistent  with  a  value  of  the  maximum  permissible 
momentum  divergence  loss,  ®div  =  0.84,  or  a  =  34-degrees.  Similar  conclusions  may  be  extracted  for  the  lower 
anode  voltage  operations. 
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propellant  utilization  efficiency  *  voltage  utilization  efficiency  =  <p(l-(?)  =  ®di»(l-(5)afi 

Fig  6.  Performance  of  Cylindrical  Hall  Thruster  (CHT)  optimum  geometry:  L  =  30  mm,  A  =  10  cm'  calculated  from  data  of  Raitses, 
Ashkenazy,  and  Guelman.2 


SQRT[Odiv(1-p)]  SQRT[®di.(1-p)] 


SQRT[<f>div(1-p)] 


Fig.  7.  Parameter  Consistency  diagrams  for  CHT  optimum  geometry:  L  =  30  mm,  A  =  10  cnr.  The  diagrams  show  values  of  fj, 
[<J>div(l-P)]1/2,  and  Q  that  are  consistent  with  each  other  at  each  of  the  six  experimental  flow  rates  and  each  of  the  three  anode  voltages: 
Va  =  300  V  (left),  250  V  (center),  and  200  V  (right). 
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VIII.  Applications  to  Hall  Thrusters  -  Comparison  of  High  and  Low  Power  Thrusters 

Figure  8  shows  a  summary  of  experimental  performance  data  (thrust  vs  specific  impulse  with  lines  of  constant 
anode  power,  and  lines  of  constant  flow  rate)  for  Flail  thrusters  ranging  from  ~  200  Watts  to  ~  50  kilowatts  and 
specific  impulse  ranging  from  ~  1000  seconds  to  5000  seconds.  The  Figure  includes  data  for  the  optimum  CFIT 
geometry  operated  at  three  anode  voltages  (Va  =  200,  250,  and  300  V)2.  Data  from  the  SP306  and  the  SP509, 
operated  between  Va  =  200  and  400  V  and  the  BFIT-200-X2B7  operated  at  Va  =  300  V  are  also  shown.  Also 
included  are  data  from  higher  power  thrusters,  P54,  173Mvl4,  173Mv24,  173MvlA24,  as  well  as  data  from  173Mvl 
with  Krypton3,  data  from  the  SP400  with  krypton  and  xenon13.  Optimum  operating  points  for  three  Busek  Flail 
Thrusters  of  (600  W14,  8  kW16,  and  20  kW16)  are  shown.  Anode  voltages  for  the  SP400  with  xenon,  Va  =  200,  300, 
400,  500,  and  600  V  are  indicated  on  the  figure  to  show  the  systematic  increase  of  specific  impulse  with  the  square 
root  of  Va.  This  approximate  dependence  of  <v>  on  Va  is  apparent  in  all  Flail  thruster  data. 


Fig.  8.  Thrust  vs  specific  impulse  for  several  Hall  thrusters  showing  lines  of  constant  jet  thrust  power  (i.e., 
jet  thrust  power  =  [F2/(2  ril  )]  =  0.1 -,  1-,  10-kW),  and  lines  of  constant  propellant  flow  rate  (0.1 -,  1-,  10-,  100- 
mg/s).  Data  from  Refs  2—1 6. 

Addition  of  the  experimental  data  on  input  power  to  the  data  of  Figure  8  produced  Figure  9,  which  is  a  primary  performance 
map  showing  the  thrust  to  power  ratio  as  a  function  of  specific  impulse.  Figure  10  shows  the  Ei  vs  E2  plot,  Ej  =  <J>(1- 
P)%  =  OdivaQ(l-(3)fi2,  and  E2  =  0(  1  -r)/y  =  (l-r)/f;Q.  Fig.  10  shows  that  the  pathways  of  increasing  f  (at  constant  r,  p.  Q,  a,  and 
Oidv)  cuts  across  lines  of  constant  qthmst,  downward  and  to  the  right,  which  is  also  the  direction  of  increasing  Va  and  m  . 

Figure  10  shows  for  the  high  power  thrusters  that  higher  anode  voltage  produced  higher  r|thmst,  which  was  accompanied 
by  smaller  (l-r)/x  in  all  the  thrusters.  It  is  interesting  that  E[  =  OdivaQ(l-P)f2  >  1  is  produced  in  SP  400  with  xenon  at 
anode  voltages  Va  >  400  V.  The  average  charge  is  the  only  factor  in  E]  that  may  exceed  unity,  so  Ei  is  the  lower  limit  of 
Q  under  these  operating  conditions  of  SP  400.. 
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Fig.9.  Performance  of  several  Hall  thrusters  with  range  of  input  power  from  200  to  50000  W  showing  lines  of 
constant  thrust  efficiency  and  lines  of  constant  specific  internal  energy  input.  Data  from  Refs.  2-16 


Fig.  10.  Performance  of  Hall  thrusters  with  range  of  input  powers  from  200  to  50000  W  with  lines  of  constant 
thrust  efficiency.  Abscissa,  Ej,  is  proportional  to  the  square  of  ion  mass  fraction.  The  direction  of  pathways  of 
increasing  f;  (at  constant  r,  |3,  a,  Q,  and  ®div)  are  shown  cutting  across  lines  of  constant  r|thrust,  approximately 
tracking  the  data  in  the  direction  of  increasing  flow  and/or  anode  voltage.  Data  from  Refs.  2-16. 
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IX.  Summary 

A  rocket  power  efficiency  equation  was  derived  to  explicitly  account  for  the  effect  of  the  velocity  distribution 
function  (VDF)  of  the  propellant  jet  on  the  conversion  of  anode  electrical  energy  to  jet  kinetic  energy.  This  enabled 
a  mathematically  rigorous  distinction  to  be  made  between  thrust  efficiency  and  energy  efficiency.  In  this  approach, 
thrust  efficiency,  the  commonly  reported  figure  of  merit  for  Hall  thrusters,  is  the  product  of  three  utilization 
efficiencies:  (1)  propellant  utilization  efficiency,  (2)  voltage  utilization  efficiency,  and  (3)  current  utilization 
efficiency.  The  three  utilization  efficiencies  are  well  behaved,  being  less  than  unity  under  all  real  operating 
conditions  and  approaching  unity  as  energy  conversion  approaches  ideality.  Perfect  propellant  utilization  is 
characterized  by  100  percent  ionization  to  a  single  ionic  species  whose  thrust  vectors  are  all  directed  along  the  same 
thrust  axis.  Anode  voltage  utilization  efficiency  is  unity  when  ion  species  are  created  at  the  anode  and  accelerated 
through  the  entire  anode  potential.  Current  utilization  efficiency  is  the  fraction  of  cathode  electron  flow  utilized  in 
neutralization  of  accelerated  positive  ions.  It  can  never  be  unity  because  a  portion  of  the  electron  flow  must  be 
recycled  back  to  the  anode  to  provide  energy  to  ionize  neutral  propellant.  The  architecture  of  the  efficiency  analysis 
is  such  that  energy  efficiency  becomes  naturally  expressed  as  a  product  of  voltage  and  current  utilization 
efficiencies.  Energy  efficiency  is  rigorously  separated  from  propellant  utilization  efficiency.  Thus,  thrust  efficiency 
is  the  product  of  propellant  utilization  efficiency  and  energy  efficiency.  The  methodology  developed  was  applied  to 
analysis  of  data  from  systematic  studies  of  performance  of  low  and  high  power  Hall  thrusters. 

Numerous  Hall  thruster  performance  studies  have  been  published  in  the  open  literature.  The  basic  data  set 
(measurements  of  thrust,  propellant  mass  flow  rate,  anode  voltage,  and  anode  current)  are  frequently  coupled  with 
measurement  of  plasma  plume  properties  using  other  optical  or  electrical  diagnostics.  Data  that  enables  calculation 
of  propellant  utilization,  <D,  and  average  charge  Q  may  be  obtained  with  an  ExB  probe,  by  Laser  Induced 
Fluroescence,  or  with  use  of  mass  spectrometry.  A  Retarding  Potential  Analyzer  provides  data  for  direct  calculation 
of  ( 1  -[I),  average  voltage  utilization.  Current  utilization  efficiency  may  be  measured  with  a  Faraday  probe.  Also, 
optical  emission  spectroscopy  may  be  applied  to  leam  about  the  thermodynamic  state  of  the  propellant  jet. 

High  power  thrusters  have  significantly  higher  propellant  and  voltage  utilization  efficiencies.  Current  utilization 
efficiency  is  similar  in  the  high  and  low  power  thrusters.  At  their  optimum  operating  points  all  Hall  thrusters  appear 
to  require  ~  10%  to  30%  of  the  total  electron  flow  to  be  recycled  through  the  plasma  to  ionize  neutrals  and  to  sustain 
dissipative  plasma  processes.  This  recycling  of  electrons  ultimately  leads  to  joule  heating  losses.  It  appears  that 
overall  thrust  efficiencies  in  the  neighborhood  of  ~  71  to  75%  that  have  been  achieved  in  higher  power  thrusters  may 
be  nearing  a  practical  ultimate  limit  because  this  requires  that  the  geometric  mean  of  the  three  utilization  efficiencies 
lie  between  89  and  91%. 
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APPENDIX 


Table  1.  Cylindrical  Hall  Thruster  [Raitses,  Ashkenazy,  Guelman,  J.  Propulsion  and  Power  14  (1998)  247-253] 


cross 
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cm2 

>> 
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300 
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0.541 
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10 

300 

1.04 

0.93 
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0.275 
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0.80 
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42.0 
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0.379 

0.647 

10 

200 

0.88 

0.61 

7.0 

0.585 

0.231 

57.9 

7.98 

0.217 

1.066 

0.371 

0.623 

6 

300 

2.16 

2.06 

33.9 

0.806 

0.43 

54.8 

15.7 

0.558 

0.771 

0.692 

0.622 

6 

300 

1.84 

1.72 

28.0 

0.791 

0.412 

54.2 

15.2 

0.525 

0.785 

0.663 

0.621 

6 

300 

1.53 

1.37 

21.4 

0.757 

0.363 

52.0 

14.0 

0.443 

0.820 

0.585 

0.620 

6 

300 

1.29 

1.28 

17.8 

0.754 

0.319 

46.2 

13.8 

0.432 

0.739 

0.573 

0.557 

6 

300 

1.04 

1.01 

13.0 

0.703 

0.267 

42.7 

12.5 

0.355 

0.753 

0.504 

0.529 

6 

300 

0.88 

0.85 

10.0 

0.670 

0.222 

39.2 

11.3 

0.291 

0.762 

0.435 

0.510 

6 

250 

2.16 

2.03 

29.5 

0.773 

0.397 

58.1 

13.7 

0.508 

0.782 

0.657 

0.604 

6 

250 

1.84 

1.67 

24.7 

0.766 

0.396 

59.0 

13.4 

0.490 

0.809 

0.639 

0.619 

6 

250 

1.53 

1.39 

20.0 

0.758 

0.378 

57.7 

13.1 

0.466 

0.810 

0.615 

0.614 

6 

250 

1.29 

1.12 

15.7 

0.743 

0.34 

56.0 

12.1 

0.401 

0.848 

0.540 

0.630 

6 

250 

1.04 

0.94 

12.0 

0.705 

0.294 

50.9 

11.5 

0.363 

0.811 

0.514 

0.572 

6 

250 

0.88 

0.74 

9.0 

0.660 

0.249 

48.7 

10.2 

0.284 

0.876 

0.431 

0.578 

6 

200 

2.16 

2.03 

25.0 

0.753 

0.357 

61.6 

11.6 

0.457 

0.781 

0.607 

0.588 

6 

200 

1.84 

1.68 

20.6 

0.744 

0.345 

61.5 

11.2 

0.428 

0.806 

0.575 

0.600 

6 

200 

1.53 

1.36 

16.3 

0.739 

0.319 

59.9 

10.6 

0.386 

0.827 

0.522 

0.611 

6 

200 

1.29 

1.05 

12.6 

0.713 

0.291 

59.7 

9.75 

0.323 

0.900 

0.453 

0.642 

6 

200 

1.04 

0.87 

9.6 

0.666 

0.256 

55.3 

9.26 

0.292 

0.878 

0.438 

0.585 

6 

200 

0.88 

0.70 

7.3 

0.639 

0.218 

52.3 

8.34 

0.236 

0.922 

0.370 

0.589 

4.5 

300 

2.16 

2.14 

32.4 

0.827 

0.379 

50.5 

15.0 

0.512 

0.740 

0.619 

0.612 

4.5 

300 

1.84 

1.75 

26.6 

0.791 

0.366 

50.6 

14.5 

0.474 

0.772 

0.599 

0.611 

4.5 

300 

1.53 

1.40 

21.4 

0.779 

0.356 

51.0 

14.0 

0.443 

0.804 

0.568 

0.626 

4.5 

300 

1.29 

1.23 

17.8 

0.765 

0.333 

48.3 

13.8 

0.431 

0.773 

0.563 

0.592 

4.5 

300 

1.04 

1.05 

13.5 

0.742 

0.279 

42.8 

13.0 

0.385 

0.725 

0.518 

0.538 

4.5 

300 

0.88 

0.89 

10.2 

0.722 

0.222 

38.2 

11.6 

0.306 

0.725 

0.424 

0.523 

4.5 

250 

2.16 

2.10 

29.5 

0.800 

0.383 

56.1 

13.7 

0.508 

0.754 

0.635 

0.603 

4.5 

250 

1.84 

1.67 

23.7 

0.774 

0.365 

56.8 

12.9 

0.450 

0.811 

0.581 

0.628 

4.5 

250 

1.53 

1.33 

19.1 

0.748 

0.358 

57.4 

12.5 

0.423 

0.845 

0.566 

0.632 

4.5 

250 

1.29 

1.12 

15.5 

0.733 

0.33 

55.0 

12.0 

0.391 

0.843 

0.534 

0.618 

4.5 

250 

1.04 

0.91 

12.0 

0.699 

0.306 

53.0 

11.5 

0.363 

0.844 

0.519 

0.590 

4.5 

250 

0.88 

0.73 

8.5 

0.632 

0.226 

46.7 

9.68 

0.255 

0.886 

0.403 

0.560 

4.5 

200 

2.16 

2.12 

24.7 

0.780 

0.333 

58.3 

11.4 

0.444 

0.750 

0.569 

0.585 

4.5 

200 

1.84 

1.66 

20.0 

0.753 

0.328 

60.2 

10.9 

0.404 

0.812 

0.536 

0.612 

4.5 

200 

1.53 

1.32 

16.1 

0.722 

0.319 

60.7 

10.5 

0.375 

0.850 

0.520 

0.614 

4.5 

200 

1.29 

1.06 

12.8 

0.708 

0.296 

59.9 

9.88 

0.332 

0.891 

0.469 

0.630 

4.5 

200 

1.04 

0.81 

9.5 

0.654 

0.268 

58.6 

9.15 

0.285 

0.941 

0.435 

0.615 

4.5 

200 

0.88 

0.61 

7.2 

0.617 

0.241 

58.9 

8.19 

0.228 

1.057 

0.370 

0.652 
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L 

cm 

<< 

m' 

mg/s 

Variable  Length  (Cross-section  =  10  crn^) 
la  F  TIP  <v> 

A  mN  X  r|  thrust  mN/kW  km/s 

Ei 

®(i-f3k 

e2 

(1-r)/x 

®(l-3) 

(1-r) 

30 

300 

2.16 

2.02 

36.0 

0.881 

0.496 

59.5 

16.7 

0.631 

0.786 

0.716 

0.692 

30 

300 

1.84 

1.67 

29.1 

0.841 

0.46 

58.1 

15.8 

0.569 

0.809 

0.676 

0.680 

30 

300 

1.53 

1.34 

21.6 

0.776 

0.379 

53.8 

14.1 

0.450 

0.842 

0.580 

0.653 

30 

300 

1.29 

1.18 

17.2 

0.745 

0.324 

48.7 

13.3 

0.402 

0.806 

0.540 

0.600 

30 

300 

1.04 

0.92 

12.6 

0.695 

0.278 

45.7 

12.2 

0.335 

0.829 

0.482 

0.576 

30 

300 

0.88 

0.80 

10.0 

0.672 

0.239 

42.0 

11.4 

0.294 

0.813 

0.437 

0.546 

30 

250 

2.16 

1.92 

31.1 

0.820 

0.466 

64.7 

14.4 

0.564 

0.826 

0.688 

0.678 

30 

250 

1.84 

1.65 

26.1 

0.803 

0.448 

63.3 

14.2 

0.546 

0.821 

0.680 

0.659 

30 

250 

1.53 

1.35 

18.9 

0.748 

0.345 

56.0 

12.3 

0.414 

0.834 

0.553 

0.624 

30 

250 

1.29 

1.11 

15.0 

0.702 

0.313 

53.9 

11.6 

0.367 

0.853 

0.523 

0.599 

30 

250 

1.04 

0.82 

11.1 

0.655 

0.289 

54.0 

10.7 

0.312 

0.928 

0.476 

0.608 

30 

250 

0.88 

0.66 

8.6 

0.618 

0.257 

52.5 

9.79 

0.261 

0.986 

0.422 

0.609 

30 

200 

2.16 

2.02 

26.1 

0.788 

0.388 

64.3 

12.1 

0.495 

0.784 

0.628 

0.618 

30 

200 

1.84 

1.66 

21.2 

0.769 

0.37 

64.1 

11.5 

0.453 

0.816 

0.589 

0.628 

30 

200 

1.53 

1.30 

15.6 

0.715 

0.305 

59.8 

10.2 

0.354 

0.862 

0.495 

0.616 

30 

200 

1.29 

0.98 

11.9 

0.682 

0.279 

60.5 

9.23 

0.290 

0.963 

0.425 

0.657 

30 

200 

1.04 

0.73 

8.7 

0.615 

0.248 

59.6 

8.33 

0.236 

1.052 

0.383 

0.647 

30 

200 

0.88 

0.60 

7.0 

0.576 

0.233 

58.2 

8.00 

0.218 

1.070 

0.378 

0.616 

40 

300 

2.16 

1.99 

34.5 

0.837 

0.462 

57.9 

16.0 

0.578 

0.799 

0.691 

0.669 

40 

300 

1.84 

1.73 

29.1 

0.827 

0.445 

56.2 

15.8 

0.569 

0.782 

0.688 

0.647 

40 

300 

1.53 

1.45 

22.9 

0.820 

0.392 

52.4 

15.0 

0.507 

0.773 

0.619 

0.634 

40 

300 

1.29 

1.21 

18.7 

0.799 

0.372 

51.3 

14.5 

0.476 

0.781 

0.596 

0.624 

40 

300 

1.04 

0.93 

14.3 

0.783 

0.353 

51.2 

13.8 

0.432 

0.818 

0.551 

0.640 

40 

300 

0.88 

0.78 

10.5 

0.700 

0.269 

44.9 

12.0 

0.325 

0.827 

0.464 

0.579 

40 

250 

2.16 

1.86 

29.9 

0.799 

0.446 

64.3 

13.9 

0.523 

0.852 

0.655 

0.681 

40 

250 

1.84 

1.64 

25.3 

0.795 

0.425 

61.7 

13.8 

0.516 

0.824 

0.649 

0.655 

40 

250 

1.53 

1.33 

19.9 

0.785 

0.389 

59.8 

13.0 

0.461 

0.844 

0.587 

0.663 

40 

250 

1.29 

1.14 

16.7 

0.755 

0.377 

58.3 

12.9 

0.455 

0.829 

0.602 

0.626 

40 

250 

1.04 

0.85 

12.7 

0.720 

0.364 

59.5 

12.2 

0.407 

0.894 

0.565 

0.644 

40 

250 

0.88 

0.68 

9.2 

0.665 

0.283 

53.9 

10.5 

0.300 

0.944 

0.451 

0.628 

40 

200 

2.16 

1.94 

25.0 

0.769 

0.374 

64.5 

11.6 

0.457 

0.818 

0.594 

0.629 

40 

200 

1.84 

1.65 

20.7 

0.769 

0.352 

62.7 

11.2 

0.429 

0.821 

0.558 

0.631 

40 

200 

1.53 

1.34 

16.6 

0.758 

0.335 

61.8 

10.8 

0.400 

0.837 

0.528 

0.635 

40 

200 

1.29 

1.02 

13.2 

0.734 

0.333 

65.0 

10.2 

0.357 

0.932 

0.487 

0.684 

40 

200 

1.04 

0.74 

10.0 

0.680 

0.325 

67.5 

9.63 

0.315 

1.030 

0.464 

0.701 

40 

200 

0.88 

0.66 

7.9 

0.620 

0.265 

59.4 

8.92 

0.271 

0.978 

0.437 

0.606 

20 

300 

2.16 

2.02 

32.9 

0.815 

0.415 

54.4 

15.2 

0.527 

0.787 

0.647 

0.641 

20 

300 

1.84 

1.61 

27.0 

0.792 

0.41 

55.9 

14.7 

0.487 

0.841 

0.615 

0.666 

20 

300 

1.53 

1.35 

21.1 

0.765 

0.358 

52.0 

13.8 

0.431 

0.832 

0.563 

0.636 

20 

300 

1.29 

1.16 

16.7 

0.712 

0.308 

47.7 

12.9 

0.378 

0.814 

0.531 

0.580 

20 

300 

1.04 

0.89 

11.5 

0.652 

0.236 

42.8 

11.0 

0.276 

0.855 

0.423 

0.557 

20 

300 

0.88 

0.71 

9.0 

0.615 

0.217 

42.4 

10.2 

0.238 

0.913 

0.387 

0.561 

20 

250 

2.16 

1.96 

29.9 

0.785 

0.424 

61.2 

13.9 

0.523 

0.810 

0.667 

0.636 

20 

250 

1.84 

1.58 

24.5 

0.759 

0.414 

62.1 

13.3 

0.483 

0.857 

0.637 

0.650 

20 

250 

1.53 

1.29 

18.2 

0.725 

0.334 

56.2 

11.9 

0.384 

0.869 

0.530 

0.630 

20 

250 

1.29 

1.08 

13.8 

0.665 

0.275 

51.3 

10.7 

0.313 

0.880 

0.470 

0.585 

20 

250 

1.04 

0.80 

10.0 

0.607 

0.241 

50.2 

9.60 

0.251 

0.961 

0.413 

0.583 

20 

250 

0.88 

0.64 

8.0 

0.570 

0.228 

50.3 

9.07 

0.224 

1.018 

0.393 

0.580 

20 

200 

2.16 

1.98 

24.0 

0.750 

0.338 

60.8 

11.1 

0.421 

0.803 

0.561 

0.603 

20 

200 

1.84 

1.57 

19.5 

0.710 

0.33 

62.2 

10.6 

0.383 

0.862 

0.539 

0.612 

20 

200 

1.53 

1.26 

15.0 

0.678 

0.293 

59.8 

9.81 

0.327 

0.896 

0.483 

0.607 

20 

200 

1.29 

1.00 

11.6 

0.661 

0.26 

58.0 

8.96 

0.273 

0.951 

0.413 

0.629 

20 

200 

1.04 

0.71 

8.0 

0.575 

0.214 

55.8 

7.67 

0.200 

1.070 

0.348 

0.615 

20 

200 

0.88 

0.58 

6.1 

0.540 

0.185 

53.1 

6.97 

0.165 

1.119 

0.306 

0.604 
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